Unconventional properties of superconducting cuprates 
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We present an explanation of the unusual peak/dip/hump features observed in photoemission ex- 
periments on Bi2212 at T < T c . We argue that these features arise from the interaction of the 
fcrmionic quasi-particles with overdamped spin fluctuations. We show that the strong spin-fermion 
interaction combined with the feedback effect on the spin damping due to superconductivity yields 
a Fermi-liquid form of the fermionic spectral function for uj < 2A where A is the maximum value 
of the superconducting gap, and a non-Fermi-liquid form for uj > 2A. In the Fermi-liquid regime, 
the spectral function A{\s.f,uj) displays a quasiparticle peak at u) — A; in the non-Fermi-liquid 
regime it possesses a broad maximum (hump) at uj 3> A. In between the two regimes, the spectral 
function has a dip at uj ~ 2A. We argue that our theory also explains the tunneling data for the 
superconducting density of states. 



PACS numbers:71.10.Ca,74.20.Fg,74.25.-q 

In recent years, the bulk of studies of cuprate super- 
conductors was focused on their unusual normal state 
properties. Less attention was paid to the behavior of 
cuprates in the superconducting state. It was generally 
believed that the superconducting behavior, even in un- 
derdoped cuprates, is rather conventional in the sense 
that most experiments can be explained in the frame- 
work of the BCS-type theory for a d— wave supercon- 
ductor. Recently, however, this belief has been chal- 
lenged by photoemission experiments on Bi2212 mate- 
rials . These experiments demonstrated that even in 
slightly overdoped cuprates, the spectral function A(k, uj) 
at T -C T c and in the momentum region near (0, 7r) where 
the d— wave gap is at maximum is qualitatively different 
from the one expected for a conventional superconduc- 
tor. Specifically, in the conventional case, A(k, u>) pos- 
sesses a single sharp peak at uj = \J A 2 . + e 2 . where Ak is 
the superconducting gap and £k is the fermionic disper- 
sion. The photoemission data for Bi2212 do show a sharp 
quasiparticle peak near the Fermi surface, but they also 
reveal two extra features in A(k, u>): a dip at frequen- 
cies right above the peak and a broad maximum (hump) 
at somewhat larger frequencies. Moreover, as one moves 
away from the Fermi surface, the sharp peak looses its 
intensity but does not disperse, while the position of the 
hump varies with k and gradually recovers the normal 
state dispersion. 

There have been several phenomenological conjectures 
in the recent literature that the sharp peak measured in 
photoemission below T c is related to the one seen in neu- 
tron scattering data in YBCO and is due to the appear- 
ance of a dispcrsionlcss mode of unknown origin below 
T c In the present communication, we present 

an alternative explanation of the photoemission data. 
We argue that the unusual superconducting properties 
of cuprates can be explained by a strong interaction be- 
tween electrons and overdamped spin fluctuations ||. 
Specifically, we show that the peak/dip/hump features 
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in the spectral function emerge due to a combination of 
two effects: (i) an almost complete destruction of the 
Fermi-liquid behavior which eliminates the quasiparticle 
peak in the normal state and gives rise to a hump in the 
spectral function at higher frequencies, and (ii) a reduc- 
tion of the spin damping at small frequencies in the su- 
perconducting state which, as a feedback effect, restores 
Fermi-liquid behavior of the spectral function in the fre- 
quency range uj < 2A. As a result, the spectral function 
near the Fermi surface possesses a quasiparticle peak at 
uj ~ A, a dip at uj » 2A where the spectral function 
experiences a crossover to a non- Fermi liquid behavior, 
and a hump at a higher frequency. As k moves away from 
the Fermi surface, the hump disperses with k while the 
quasiparticle peak only looses its intensity as it cannot 
move farther in frequency than 2A. This behavior fully 
agrees with the photoemission results jl J~ 

We first briefly review the results for tl 
tion in the normal state and then discuss our calculations 
in the superconducting state. One of us has recently 
shown H that the fermionic self-energy due to the spin- 
fermion interaction is almost independent of momentum 
and has the form g 2 T,(uj) — g 2 2uj/(l + \/l ~ i\uj\/uj s f) 
where g 2 cx £ is a dimensionless coupling (£ is the 
magnetic correlation length) and uj s f oc Q? £) is the 
typical relaxation frequency of overdamped spin fluc- 
tuations. It was further argued that in cuprates g 2 
is large even for overdoped materials so that the self- 
energy overshadows the bare frequency dependence in 
the fermionic propagator, i.e., G„(k,w) w Z/(Y>(uj) — £k) 
where Z — g~ 2 , £k = Ze^, and ek is the bare fermionic 
dispersion. For uj < uJ s f, S(w) has the Fermi liquid 
form Tt(uj) uj + iuj\uj\/ (Auj s f) and accordingly, the spec- 
tral function A(k, uj) has a conventional peak at uj = 
though with a reduced residue Z. For uj > uj s f, however, 
the system crosses over to a region which is in the basin of 
attraction of the quantum critical point, £ = oo. In this 
region, g 2 ^(uj) rj ^4e l7r / 4 |o;| 1 / 2 sgn(w) where A = 2g 2 u) 1 Jj 1 
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FIG. 1. The calculated quasiparticle spectral function in 
the normal state (a) and in the superconducting state along 
M -T ((0,tt) - (0,0)) (b) and M - Y ((0,tt) - (tt.tt)) (c). 
The results are presented for 6 = 5. For the M point, we used 
e = A. 

is independent of £. As a result, instead of a sharp 
quasiparticle peak, the spectral function possesses only 
a broad maximum at co = e^/(4w s f) ^> e k - Various 
NMR experiments have shown that co s f is anomalously 
small even if the correlation length is comparable to the 
interatomic spacing (co s f < 15 meV even at optimal dop- 
ing |]l0|) Since this lo s j is smaller than the resolution of 
the photoemission experiments, the experimentally mea- 
sured A(k, co) only displays a broad maximum even for 
k k, k F (see Fig. |l|a,). 

Consider now the same system at T « T c . We ar- 
gue that there are two key effects associated with the 
superconducting state. First, the quasiparticle Green's 
function is modified due to the fermionic pairing. In the 
BCS approximation we have Gj c 1 (k,w) = G,7 1 (k, co) + 
A k /Z) 2 G„(-k, -co) where (A k /Z) 2 is the strength of 
the d-wave pairing susceptibility. In a Fermi-gas (Z = 1) 
Afe is the pairing gap. In our case, however, Z < 1, 
G n (k, oj) = Z/(Y>(co) — efe), and we obtain 

G sc (k,o,) = Z s2M _ (A , +4) . (1) 

For co < co s t, this form of G sc (k, to) again yields a conven- 
tional quasiparticle peak at co = yj A 2 . + e 2 ,. For co > co s f , 
however, the spectral function takes the form 

Arir \ /n M + E k + € kv / \co\/{2co sf )sgn{co) 

4 ( k ' U ) K VM ^2 + E 2 ( 2 ) 

where Ek = (A 2 , + e 2 ,)/(4tj s /), and thus shows the same 
features as the one in the normal state: it possesses a 
broad maximum at co = Ek, but no quasiparticle peak. 

Clearly, the possibility to observe a quasiparticle peak 
in cuprates near k = (0, it) where A k is near its maximum 



value A depends on the ratio b — A/co s f. If this ratio is 
small, then at frequencies comparable to A, the system 
is in the Fermi liquid regime and A(k, co) displays a peak, 
while if b > 1, then at lo ~ A, the system is already in 
the non-Fermi liquid regime where AQc, co) only exhibits 
a broad maximum. An earlier computation of b in the 
spin-fluctuation approach || has shown that b increases 
with decreasing doping and becomes larger than one al- 
ready for slightly overdoped cuprates. At this stage, we 
therefore only obtain a hump in A(k, co). 

However, there is a second effect related to supercon- 
ductivity which gives rise to a quasiparticle peak in the 
spectral function even if 6 ^> 1. Indeed, in the above 
analysis we have so far assumed that co s f is independent 
of frequency. This, however, is not true in the super- 
conducting state as the opening of the superconducting 
gap reduces the spin damping at low frequencies and 
hence increases co s f in the same frequency range. From 
this perspective, the parameter b we introduced before 
is identical to boo = A/co s f(co = oo) where co s f(co = oo) 
is the value of co s f in the normal state, while the exis- 
tence of a quasiparticle peak is actually determined by 
b u = A/co s f(co) at co ~ A. 

The spin relaxation frequency co s f is inversely propor- 
tional to the damping of a spin fluctuation at Q = (it, n) 
due to its decay into a particle-hole pair. To obtain 
co s f(co) we hence have to evaluate the imaginary part of 
the particle-hole bubble. The normal state calculations 
have been reported previously p] and yielded a frequency 
independent co s f = (3/16)Zvf£~ 1 ( Vp is the Fermi ve- 
locity). It is essential that this result is independent of 
S and is therefore the same as for free fermions. This 
universality is a general consequence of the fact that the 
normal state self energy depends only on frequency ||. 
Vertex corrections do modify co s f, but these corrections 
are very small numerically and can be safely neglected || . 

We now turn to the superconducting state. Here we 
have to (i) combine two fermionic bubbles made of nor- 
mal and anomalous Green's functions and (ii) reevalu- 
ate the fermionic self-energy S(k, co) using superconduct- 
ing Green's functions for intermediate fermions. Simple 
estimates show that for finite A, co s f indeed acquires 
a frequency dependence and recovers the normal state 
value only at lo 3> A. Moreover, unlike normal state 
G^ x (k, co), superconducting G ; T c 1 (k, co) cannot be written 
as the normal state Fermi-gas result plus a momentum- 
independent self-energy. As a result, co s f(co) depends on 
the self-energy of the intermediate fermions which by 
itself depends on co s f. Consequently, the equation for 
co s f(co) becomes an integral equation. Below we solve 
for co s f(co) assuming that T,(k,co) in the superconducting 
state has the same form as in the normal state albeit with 
a frequency dependent co s f(co). We checked by explicitly 
evaluating the lowest-order self-energy diagram that the 
terms we omitted by using the normal state rather than 
the superconducting Green's function change the results 
by less than 5% for all frequencies. 

To obtain the equation for b u , we integrated over the 
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fermionic momentum in the bubbles and assembling nor- 
mal and anomalous contributions to the spin damping. 
This yields 



- Re 

UJ 



A 2 - E(ft+)£(ft_) + D{n + )D(n_) 



D(n + )D(n_) 



(3) 



where, we remind, — A/u) s f(u>), b = b^ = A/ uj s f(oo). 
Also 0± = n±w/2, £(fi) = 2X1/(1 + y/l~i\fl\/uj s f(n)), 

and = V s2 ( n ) - A2 - 

It is instructive to consider the limits of small and 
large b. For d « 1, the fermionic damping is negligi- 
ble, £(Q) « fi, and 6 W has the same functional form as 
the spin damping in the superconducting Fermi gas: it is 
almost zero for uj < 2A, jumps to 6 W = bir/2 at uj > 2A 
and gradually approaches b with increasing frequency. A 
more careful look into Eq.(||) shows that the jump in 6 W 
at uj — 2A is a consequence of a singularity in the in- 
tegrand at |fi±| = A. For finite fermionic damping, the 
singularity is washed out but as long as the damping at 
|f2±| = A is small (which is clearly the case for small 
b), b^ sharply drops below u> = 2 A. Consider next what 
happens when b S> 1. For uj 3> A, uj/uj s f(uj) 3> 1 and 

hence E(fi ± ) oc e^Vl^il- Using this form of £(fi±), 
one can easily verify that the downturn renormalization 
of b u starts already at ui ~ Ab jP A. At smaller frequen- 
cies, we expand the integrand in ([}]) in powers of uj and 
obtain b u ~ oj/A. For ~ A, we then have b u ~ 1 mde- 
pendently on how large & is. This in turn implies that at 
these frequencies, the system crosses over into the Fermi 
liquid regime, and the behavior of b^ at uj < A is quali- 
tatively the same as for small b. Clearly then, 6 W should 
drop at around uj = 2A no matter how large b is. Our 
numerical solution of (^) which we present in Fig. || con- 
firms our analytical result. We see that while the high 
frequency part of b^ evolves with increasing b from fre- 
quency independent to a linear in frequency behavior, a 
sharp drop in b u survives even for large b. 

The strong reduction of b^ at low frequencies yields a 
crossover in the system behavior from a non-Fermi liq- 
uid, strong coupling behavior at uj > 2A to a conven- 
tional Fermi-liquid behavior at uj < 2A. In the latter 
case, the spectral function A(k.F,uj) indeed possesses a 
sharp quasiparticle peak at uj = A and rapidly decreases 
at higher frequencies. Combining these two limits, we 
find three distinct features in A(\tp,uj): a quasiparti- 
cle peak at uj ~ A, a dip at uj ~ 2A and a hump at 
uj = A6/4 which at 6 ^ 1 is parametrically larger than 
the peak frequency (see Fig. [l]). These results are fully 
consistent with the photoemission data [QJ|]. Further- 
more, analyzing the form of the spectral function Eq.(jl]) 
for various momenta, we find that the quasiparticle peak 
does not disperse with k as the region of Fermi liquid 
behavior does not extend farther than 2A away from the 
Fermi surface. Instead, the peak gradually decreases in 
the amplitude as one moves away from kp. This is clearly 
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FIG. 2. The solution of the integral equation for the fre- 
quency dependent spin-fluctuation frequency w s f(u>) in the 
superconducting state. In the normal state, ui s f = oi s /(oo). 
For all b = A/u s /(oo), ujJ}(uj) sharply drops at around 2A. 

seen in Fig. [l]b,c. In contrast, the position of the hump 
follows uj tx (A 2 , + e|). Near (0, 7r), where ek is small, 
the dispersion is weak, whereas further away from the 
Fermi surface it disperses with k and gradually recovers 
the normal state dispersion. In Fig |3[ we plotted the fre- 
quency position of the quasiparticle peak and hump in 
the superconducting state vs. the normal state position 
of the hump. This dependence on momentum is also fully 
consistent with the photoemission data |jj . 

Next we compute the density of states N(uj) in the 
superconducting state assuming that the dominant con- 
tribution to N(uj) comes from momenta near (0, 7r). In- 
tegrating A(k, uj) over ek, we obtain 



N(uj) cx Re 



E(w) 



V^ 2 (oj) - A 2 



(4) 



The plots of N(u>) for various b are presented in Fig. |J. 
We see that for all values of b, N(u>) possesses a peak at 
uj = aA where a sa 1 for small b and gradually increases 
with increasing b. For larger &, N(uj) displays a dip at 
frequencies slightly larger than 2A; the amplitude of the 
dip increases with b. Above the dip, N(uj) increases as 
y/uj and eventually saturates. These results are in full 
agreement with the tunneling data in Ref. except for 
the experimentally observed anisotropy between N(uj) for 
positive and negative frequencies for which we do not 
have an explanation. 

Furthermore, we find that for large b (i.e., for under- 
doped cuprates), the region near (0, it) yields a substan- 
tial contribution to N(uj) even for uj < A, which possi- 
bly overshadows the contribution from the nodes of Ak- 
Specifically, we found from Eq.(|]) that N(uj) oc uj 3 for 
very small frequencies which crosses over to N(uj) oc uj 
around uj ~ A/2. This behavior is very similar to the one 
observed in underdoped cuprates. On the other hand, for 
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FIG. 3. The frequency position of the quasiparticle peak 
and the hump in the superconducting state versus position of 
the hump in the normal state obtained from Fig. |tj. 
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FIG. 4. The tunneling density of states in the supercon- 
ducting state. Observe the development of a dip with in- 
creasing b. 



small b (i.e., for strongly overdoped cuprates), the contri- 
bution from the nodes is indeed relevant and our results 
underestimate N(uj) for \u\ < A. 

Finally, we discuss the value of the gap and its variation 
with doping. Consider first the location of the quasipar- 
ticle peak in A(k,u>). Our theory predicts that at some 
distance away from the Fermi surface, it should be lo- 
cated at oj ss 2A. Applying this result to nearly optimally 
doped Bi2212 materials studied in Ref. we obtain 

A ~ 25 — 30 meV. Almost the same result is obtained 
by extracting 2A from the onset of the dip in measured 
N(uj) 0. Furthermore, for T c = 83K underdoped mate- 
rial, our analysis of the tunneling data yields almost the 
same value of A as at optimal doping. We therefore ar- 
gue that A almost saturates around optimal doping, and 
the observed increase of the peak frequency in N(co) with 
decreasing doping is mostly due to strong coupling effects 
which shift the peaks towards higher frequencies. Notice 
also that a comparison of A with uj s f(oo) extracted from 
NMR data in the normal state confirms our assertion 
that b > 1 already at optimal doping. 

To summarize, in this paper we present the explana- 
tion of the unusual peak/dip/hump features observed in 
photoemission experiments on Bi2212 at T <C T c . We 
argue that these features are explained by the interac- 
tion of fermionic quasiparticles with overdamped spin 
fluctuations. We show that the strong spin-fermion in- 
teraction combined with the feedback effect on the spin 
damping due to fermionic pairing yield a Fermi-liquid 
form of the fermionic spectral function at w < 2A where 
A is the maximum value of the superconducting gap, 
and a non-Fermi-liquid form for to > 2A. In the Fermi- 
liquid regime, the spectral function near kp displays the 
quasiparticle peak at lu rs A; in the non- Fermi-liquid 
regime above 2A it possesses a broad maximum (hump) 
at u> ^> A. In between the two regimes, the spectral func- 



tion has a dip at uo ~ 2A. We argue that our theory also 
explains the tunneling data for the superconducting den- 
sity of states. We predict that the superconducting gap 
saturates around optimal doping, and that the observed 
increase of the peak frequency in the tunneling density 
of states with decreasing doping is chiefly due to strong 
coupling effects. 
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